#include <iostream>

using namespace std;
const int N=1e7+10;
const int MOD=1e9+7;

int fac[N];
int invfac[N];
void init()
{
	fac[0]=1;
	for(int i=1;i<=N;i++)
	{
		fac[i]=(fac[i-1]*i)%MOD;
	}


}
int main()
{
	int q;
	cin>>q;
	init();
	while(q--)
	{
		int n,m;
		cin>>n>>m;
		cout<<(fac[n]-fac[n-m])/fac[m]<<endl;
	}	
	return 0;
}
